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1. Introduction 

Recently, we adopted the dynamic theory of uniaxial liquid crystals (Ericksen- 
Leslie's theory) to a two-dimensional thin film so as to make the theory applicable 
for the study of hydrodynamic behaviour of lipid bilayer membranes [1] . The mem- 
brane was considered to be isolated from spatial surroundings, which is evidently 
impossible to realize in experiments. In aqueous solutions lipid membranes tend to 
close and form vesicles, due to the spontaneous curvature, maintaining a difference 
of pressure between the interior and exterior liquids. Under the action of chemical 
or electrochemical fields, solutes of the adjacent solutions may migrate across the 
membrane, inducing oscillations of the membrane potential [2-4], whereas under 
t^j- ! the action of physical fields, shape or structure of the membrane may change and 
may keep varying for a long duration [5-8]. Such mass or momentum exchange 
; between the membrane and its surroundings is, indeed, essential for the complex 
dynamic processes. The purpose of this paper is to connect in theory the film of 
liquid crystals with the bulk phases of normal fluids. 

The theory will be extended to the membrane-solution system in the following 
way: one omits the thickness of the membranes, as in the foregoing paper [1], 
and thus is able to regard the two bulk phases as if they were in direct contact 
with each other. It follows that the "dividing surface" of the two liquids coincides 
with the "middle surface" of the membrane. In any irreversible process, on the 
one hand, kinematic discontinuity occurs at the "dividing surface" of the bulks, as 
it is expressed by the condition of compatibility (Kotchine's theorem [9]); on the 
other hand, instantaneous balance of mass, momentum and energy is required at 
the "middle surface" of the membrane, as it is presented in the context of isolated 
membranes [1]. It stands then to reason that a superposition of these individual 
laws gives the description of the desired dynamical system. It is done in the second 
section. A group of dynamic equations for an open (and closed) membrane system 
is thus established. 

Here we explicitly note that the lipid crystal membrane is a two-dimensional 
fluid, which responds to lateral forces with the surface tension (or compressibility) 
and to transverse forces with the elasticity, having no elastic response to shear 
stretch. This is in contrast to the solid film model for red cell membranes suggested 
by Evans [10]. In the gel state, the lateral fluidity vanishes and the shape of the 
membrane varies in the bending mode (if no molecule tilting). 

Recent experiments give support to the adoption of liquid crystal theory to 
the study of biomembranes. It was found that membrane deformations, at least 
those of small dimension such as shape fluctuations, are dominated by the curva- 
ture elasticity and are essentially "shear- free" . The same conclusion was obtained 
from the observations on both artificial lipid vesicles [11] and on erythrocytes [12]. 
Another argument is given by the measurements of the bending modulus (k'n). 
It was reported that the bending modulus of lipid membranes is very close, in 
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magnitude, to that of human erythrocyte membranes: for phospholipid vesicles 
fc; i = l~4x 10" 12 erg [5,13,14]f and for red blood cells k' lx = 1 ~ 7 x 10" 12 erg 
[12,15,16]|. 

We expect that the theory presented in this paper will be extensively used in 
the study of dynamic processes of lipid vesicles, as well as biological cells, suspended 
in aqueous solutions. As a special application, one gives in section 3 the hydrostatic 
equations of the membrane in the field of fluid pressure and suggests to interpret the 
complex cell shapes by means of the analysis of stability, since biological systems 
are in general far from thermodynamic equilibrium. 



2. Connection with liquid surroundings 

The bilayer lipid membrane is considered to be adjacent to a stokesian fluid 
on each side. The thickness of the membrane is neglected and thus the "middle 
surface" of the membrane coincides with the "dividing surface" of the two fluids 
(see Fig. 1). The position of the two superposed geometric surfaces in an Euclidean 
space is given by the vector 

r 1 = r 1 {e\ 2 , t) (1) 

where 6 a are curvilinear coordinates on the surface and t is the time. Here and 
throughout the paper, block letters stand for vectors or tensors expressed in the 
inertial coordinates fixed in the Euclidean space and italic letters for those in the 
local bases (see reference [1]); Latin suffixes range over 1, 2, 3, while the Greek ones 
over 1, 2. 

Denoting the displacement of the surface by £, one defines the speed of the 
displacement movement by 

~£ < 2 > 

Both C and w are measured in the direction normal to the surface (see the definition 
of the displacement movement [1]). 

2.1 Mass balance 



f The data cited here are for egg lecithin. For pure phospholipid membranes, it 
was found that k'n = 4 ~ 5 x 10 -13 erg [17,18]. Adding cholesterol to the pure 
lipid membranes may elevate the magnitude to one order higher. For example, the 
modulus of dimyristoylphosphatidycholine containing 50 mol% cholesterol has the 
value k'n = 4 x 10 -12 erg [19]. This composition equals to that of erythrocyte 
membranes which contains 50 mo/% phospholipids and 50 mo/% cholesterol. 

| The data cited here correspond to the exciting wave length comparable to the 
cell dimension (10 fim). It was also reported that to short wave length excitation 
(0.1 ~ 0.3 |im), the bending modulus of red cell membranes is 1 ~ 2 x 10~ 13 erg 
[20,21]. For a discussion about the discrepancy refer to [12]. 
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In the absence of exchange with the surroundings, the membrane has a constant 
total mass. The mass conservation principle for an arbitrary area S is written 



d_ 

dt 



7cte = (3) 



s 



where 7 is the surface density of mass. 

When the membrane is permeable to certain solutes but not to the solvent, the 
total mass balance is given by 

J t JJ^ ds = ~J2jJ b'fe(n)] ds (4) 

S k S 

where jfc is the current density of the kth considered species crossing the membrane. 
The subscript (n) denotes the normal component of vectors. The square brackets 
stand for the discontinuities. For an arbitrary physical quantity 4> the jump at the 
surface is defined as 



with the subscript cone for the concave side and conv for the convex side of the 
surface. 

When osmosis of the adjacent solutions occur, the instantaneous mass change 
is additionally furnished by the absorption flux over the area. So one has 

J t fjl ds =~Y^ jf [h(n)] ds ~ j j [P ( U (n) ~w)]ds (5) 

s k s s 

where p is the mass density and u the velocity of the bulk fluids. 

However, the Reynolds transport theorem for two-dimensional fluids asserts 
that [22] 

d_ 

dt 



4>ds 



ds (6) 



where, 4> is any scalar function of position and time, a the determinant of the 
metric tensor of the curvilinear coordinate system and v the velocity of the fluidsf. 
The top dot stands for the convected derivative and the comma for the covariant 



f In this paper, u denotes the velocity of the surrounding fluids, whereas v 
denotes that of the membrane. In the local bases, the velocity of the membrane 
(without thickness) is expressed as 



v a e a + we 3 



with e3 = n, n being the normal unit vector. 
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derivative along the surface. Summation convention is followed here and throughout 
the paper. Besides, the dilation of the area element is related to the displacement 
movement by [9] 

A = - 2Hw (7) 

H being the mean curvature of the surface. So Eqs.(3-5) may be written in the 
differential form 

j + j(v« a -2Hw) =0 

(for closed membrane system) (8) 

i + l{v a a -2Hw) = - J £[jk(n)} 

k 

(for solute permeation process) (9) 

7 + 7 (v a a ~ 2Hw) =-J2 [jk(n)] ~ [p («(„) - w)] 

k 

(for solution osmosis process) (10) 

Here we have omitted the lateral diffusion of the particles moving across the mem- 
brane. 

2.2 General form of balance equations 

Suppose first that no exchange of the physical quantity occurs between the 
membrane and the bulk phases. Let <p = 7<^ s - One may write then the instantaneous 
balance in the form 

j t JJ^ s ds = JjQds- jj a dl a (11) 

S S L 

where L is the periphery of the area S, Q the production of <p pe r unit time per 
unit area and J the current density of <fi along the surface. Making use of theorem 
(6), as well as the Stokes theorem, and taking account of (7), one obtains 

4 8 + <f [7 + 7 (<« - 2Hw)] =Q~J a a (12) 

Imagine next that the two bulk phases are in direct contact with each other 
without the presence of the membrane between them. The condition of compatibil- 
ity at the dividing surface is given as (Kotchine's theorem [9]) 

[p4> v (u (n) - w)} + [j (n) ] = (13) 

where pcff = (p and j is the influx of <f> from the volumes to the surface. 

Then we consider the membrane being in contact with the surroundings. Evi- 
dently, in this case the accumulation brought by the flux from the volumes makes a 
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part of the total amount of <fi pertinent to the membrane. As the middle surface of 
the membrane and the dividing surface of the two bulks are coincident, one obtains 
naturally the desired balance equation by summing Eqs.(12) and (13) 

7 <^ + s [j + 7 ( v a ja -2Hw)] 

= Q-J%- [3(n)} ~ W («(„) - W)} (14) 

In particular, put (j) s = (p v = 1, then Eq.(14) will be reduced to Eq.(lO), if the 
source of substances is null and the lateral diffusion of components is negligible. 

Inserting (8), (9) and (10) respectively into Eq.(14) gives the balance equations 
of </>, for the closed membrane system 

4 s =Q - J% - [j [n) ] - [ P 4> v («(„) - w)] (15) 
for the solute permeation process 

4 S =Q - - [j {n) ] - [pfr {U {n) - W)] +^ S J2 [jk{n)] (16) 

k 

and for the solution osmosis process 

4 S =Q - - [j{n)\ ~ W («(n) - W)] + 

+ S J2 b'fc(n)] + S [P i U (n) ~ W )} ( 17 ) 
fc 

In the following, one considers only the closed membrane system. 

2.3 Momentum exchange 

With regard to the momentum exchange, the physical quantity <f> is interpreted 
as the velocities 

s = + W 



(j) v = u 1 (18) 

where n is the unit normal vector to the surface and t is the 3x2 hybrid tensor 
having the components 

= w° (19 » 

x 1 being the inertial coordinates. The net influx [j( n )] is interpreted as the resultant 
of the forces that the bulk phases exert on the surface 

[j(n)} = [q] g^n 1 - [pW + 2 [^e ij ] nj (20) 

where g equals to the mass density p multiplied by a unit length, g is the gravita- 
tional acceleration, p the thermodynamic pressure, e^ the rate of strain tensor and 
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H v the viscosity coefficient of the bulk fluids. The source of is interpreted as the 
gravitational attraction to the membrane 

Q = 7S 1 (21) 

The surface flux J is interpreted as the surface force inherent in the membrane 

J a = -T ap t) 4 - T a3 n [ (22) 

where T is the stress tensor of the membrane. With these interpretations, Eq.(15) 
reads 

lj t + urn 1 ) 



= ^+(T a % + T a3 n i ) j 

-[g]g^n i + \p}n i -2[^}n i (23) 

In writing (23) we have supposed the continuity of velocity at the surface 

U(n) = w (24) 
Making use of the kinematic formulae [23,24] 

4 = (W) a + 4V (25) 



, Of 

and 

h 1 ^n' + n 1 ^ + w j / k jn j n k 

= -t i a (a af3 dpw + b%v f3 ) (26) 
where a a p and b a p are respectively the first and the second magnitude of the surface, 



i 



is the second kind of Christoffel symbol, one rewrites the first term on the 



jk 

left side of (23) as 



= 7 - 2b%v f3 w - a^wdpw) t l a + 

+ 7 (w + v a d a w + b aP v a v p ) n 1 (27) 

Additionally, one expands the second term on the right side of (23) in the form (see 
Appendix) 

= (T^ {T) - b^ 3 ) % + (T«l (T) + b af) T#) n 1 (28) 



with 



T ap , , =d a T ap + —T a(3 d a a +{ & \t c 
- a ( T ) 2a [a a) 



™1 {T) =d a T a3 + ^T a3 d a a (29) 

where the subscript (r) is used to denote the tangential component of vectors or 
tensors. Then, by using (27) and (28), one splits Eq.(23) into the tangential part 

7 (i) a - 2b a p v p w - a a(3 wdpw) 

= 7 gj*X a + Tf J (T) - b%T& - 2 [^ eij ] MpaP a (30) 
and the normal part 

= 7 gjuj + Tf| (T) + bp a T a ^ - [q] gj nj + [p] - 2 [^,ey] (31) 
The rate of strain tensor for the stokesian fluids is given, as well known, by 

ey = ^Kj +Uj,i) (32) 

The stress tensor of the membranes without any chiral structure has the components 
[1] 

T af3 = - aa af3 + (k[ - 2Hk , 11 ) b af3 + 

+ (k' 22 + k' 2A ) Ka ap -(jjl-ti) 2Hwa ap - 



T a3 = _ ?* b c< v (33) 

where, the k'^s and the fe^-'s are the elastic moduli, \i and rj the dilatation vis- 
cosity and the shear viscosity, respectively, a is the surface tension, ct 3 the viscos- 
ity relevant to the intrinsic rotation of the molecules, K the gaussian curvature, 
e" a the two-parametric permutation tensor having the components e 11 = e 22 = 
0, e 12 = — e 21 = l/v^a, _E /3aA/i the fourth-order transversely isotropic tensor de- 
fined as E 13 ^ = a aX a 0fl + a afl a f3x - a a(3 a x ^. Inserting relations (32) and (33) into 
Eqs.(30) and (31), one obtains finally the explicit form of the equation of motion 
for the internal flow within the membrane 

7 (i) a - 2b a v p w - a ap wdpw) 

= 7g j t> /3a - aP a dpo + k' x aP a dp (2H) - 

- k' u (b pa + 2Ha Pa ) dp (2H) + (k' 22 + k' 2A ) a pa dpK- 
-(ji + rj) a Pa dp (2H) w - {(/i - rj) 2Ha ap + 2r]b ap ) dpw+ 
+ \E^ (v x , m/3 + X p) + ^ {2Hb%vP - Kv a ) - 

-KKj+Uy)]!!^ (34) 

and that for the displacement movement of the membrane 
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7 (w + v a d a w + b a pv a v 13 ) 

= 7g jn j - 2Ha + (k[ - 2Hk[ 1 ) {AH 2 - 2K) + (k' 22 + k' 2A ) 2HK- 
-4{(p + V )H 2 - V K}w + V (b x » - Ha x ») (v x , „ + v», x) - 

- y W v °p + v % ^} ~ W ^ + W- 

- [n v Kj +u j)i )]n 1 n J (35) 

2.4 Energy exchange 

With regard to the energy exchange, one might interpret the physical quantity 
4> in (15) as the energy of unit mass, the surface flux J as the sum of the current 
density of heat and work, and so on, to get the balance equation. But here we use 
the results obtained before to make the discussion brief. 

One has already the energy balance equation of an isolated membrane [1] 

-yQQ = Q h - J h * a + T^ a v a , p + T&dpw (36) 

with the constitutive relation of heat flux 

J h = -KL^e^e" - K||djen j e 3 (37) 

where, C\ is the heat capacity per unit mass of the membrane medium in the static 
state, © the temperature, Q h the heat source, k the heat conductivity with the 
subscripts _L or || indicating the conduction direction perpendicular or parallel to 
the normal vector and e k are the reciprocal local bases. While the membrane is in 
contact with the bulk fluids, energy current flows from the volume phases to the 
surface in the form of heat and work 

- [j(n)] = [K v Q,k] n k - w {[q] g (n) - [p] + 2 [/^eij] n'n J } (38) 

where k v is the heat conductivity of the volume phases. By summing (36) and (38), 
one obtains the energy balance equation 

7 C/6 =Q h + K± (<9je) k (g k J - n k nJ) - n a 2H (n^e) + 

+ (^-aa ap + (k[ - 2Hk' ll ) b a0 + 

+ (k' 2 2 + k' 2 d Ka ap -(jjl-ti) 2Ha ap w- 

- V 1 2wb^ - (v Xt „ + « M , a) }) v a ,p- 

~ Y^ 9 ^ + i^O] n k - 

- w\[q] gj nJ - [p] + [u. v ( Ui)j + Uj)i )] n'nU (39) 
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where K a = k\\ — k±. In writing (39), one took relations (33) and (38) into account. 



3. Hydrostatic equilibrium 

Let v a = w = u 1 = 0, Eqs.(34) and (35) give 

7g • t - V s a + k[V s (2H) + (k' 22 + k' 2A ) V S K+ 

+ k' n {(V s n) • V s (2H) - 2HV S (2H)} = (40) 

and 

7g(n) - 2Ha + k[ {AH 2 - 2K) + (k' 22 + k' 24 ) 2HK- 

- k' u 2H {AH 2 - 2K) - [ e ]g (n) + [pj = (41) 

where the differential operator is defined by 

V. = e^ (42) 

For a non-spherical vesicle, the pressure difference between the interior and the 
exterior may be estimated to be of the order of 1 ~ 10 dyn ■ cm~ 2 [25], whereas 
the gravitational attraction 7g is approximately 5 x 10~ 4 dyn ■ cm~ 2 . So the later 
is negligible as an external force. On the other hand, theoretical analyses and 
experiments show that the saddle-shape conformation, given by the term of k 22 + 
k' 2A , can not be observed in lipid- water systems except when the water content is 
extremely low [26]. Therefore Eqs.(40) and (41) can be simplified to 

V s a - k[V s (2H) - k' n {(V s n) • V s (2H) - 2HV s (2H)} = (40') 

and 

2Ha - k[ {AH 2 - 2K) + k' n 2H {AH 2 - 2K) + [g]g (n) - [p] = (41') 

Equations (40') and (41') serve as the conditions of mechanical equilibrium for 
the membrane-fluid system in the directions tangential and normal to the surface, 
respectively. Eq.(41') gives explicitly the dependence of the elasticity of the mem- 
brane on the discontinuity of the pressure field. It is worth while to note that these 
conditions are different from those of thermodynamic equilibrium [25,27,28]. For 
the latter requires further the minimality of the free energy stored in the whole elas- 
tic surface. When the elasticity vanishes and the buoyancy is negligible, Eq.(41') is 
reduced to the well known Laplace law 

2Ha = [p] (43) 

4. Concluding remarks 

Eqs.(15), (16) and (17) are given as the general form of balance equations for 
membrane-solution systems. Specifically, Eq.(15) is for the process without mass 
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exchange between the membrane and the liquid surroundings, Eq.(16) for the pro- 
cess of solute permeation and Eq.(17) for the solvent osmosis across the membrane. 
Particular uses of Eq.(15) yield the equations of motion (34) and (35), respectively 
for the internal molecular flow inside the membrane and the displacement of the 
membrane configuration, as well as the balance equation of energy Eq.(39). The 
conditions for the membrane in hydrostatic equilibrium with the adjacent fluids are 
given as a special application of Eqs.(34) and (35). 
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Appendix. Derivation of Eq. (28) 

Recalling the covariant differentiation [22] 



a/3 d6P \ j k J a P \a(3 
one writes the first term of (28) as 

(T a %) , a =d a (T a %) + { ^ } t»% + { .\ } 4#r* 

and the second term as 

(r«V) a =d a (r«V) + 1 .\ } i^r-V + 1 a « } r"V 

=n' (s„r- + { a « } r-*) + 1- (y + { / fc } 

_ ima3 , t«3 i 
— n J ,a(r) + 1 n ,a 

~ n J ,a(r) J a r /3 

Hence one obtains 

(T a/3 ^ + T a V) q 

= (^f (T) - T^b£) % + (T«l (T) + T«% a ) n 1 
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Figure captions 

Fig.l Illustration of the middle surface of the membrane and the dividing sur- 
face of the two bulk phases. 
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